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Abstract. We present the first circular polarization mea- 
surements of circumstellar H2O masers. Previously the 
magnetic field in circumstellar envelopes has been esti- 
mated using polarization observations of SiO and OH 
masers. SiO masers are probes of the high temperature 
and density regime close to the central star. OH masers 
are found at much lower densities and temperatures, gen- 
erally much further out in the circumstellar envelope. The 
detection of the circular polarization of the (615-523) rota- 
tional transition of the H2O maser could be attributed to 
Zeeman splitting due to the magnetic field in the interme- 
diate temperature and density regime. The fields inferred 
here agree well with predicted values for a combination of 
the r~^ dependence of a solar-type magnetic field, and the 
coupling of the field to the high density masing regions. 
We also discuss the unexpected narrowing of the circular 
polarization spectrum. 



1. Introduction 

High mass loss in late type stars produces a circumstellar 
envelope (CSE) in which several different maser species 
can be found. These masers, especially SiO, H2O and OH, 
are excellent tracers of the dynamics and kinematics of 
the CSEs. Polarization observations of these masers have 
revealed the strength and structure of the magnetic field 
throughout the CSE. Observations of SiO maser polariza- 
tion have shown highly ordered magnetic fields close to 
the central star, at radii of 5-10 AU where the SiO maser 
emission occurs (Kemball & Diamond, 1997). The stan- 
dard Zeeman interpretation gives magnetic field strengths 
of ~ 5-10 G. However, a non-Zeeman interpretation has 
been proposed by Wiebe & Watson (1998), which only re- 
quires fields of « 30 mG. At much lower densities and tem- 
peratures and generally much further from the star, OH 
maser observations measure fields of w 1 mG (Szymczak 
& Cohen, 1997; Masheder et al., 1999). But for the inter- 
mediate region, at distances of a few hundred AU, no in- 



formation is available. In this region the H2O maser emis- 
sion occurs. Since water is a non-paramagnetic molecule, 
determination of the magnetic field is significantly more 
difficult. The Zeeman splitting of H2O is extremely small 
for the field strengths expected (few hundred mG), only 
« 10^'' times the typical Gaussian line width of the H2O 
maser line (Ai^l ~ 20 kHz). However, Fiebig & Giisten 
(1989, hereafter FG) showed that in the presence of such 
magnetic fields the Zeeman splitting can be detected with 
high spectral resolution polarization observations. Their 
observations targeted strong interstellar H2O maser fea- 
tures. The observations presented here give the first results 
of circular polarization measurements of the H2O masers 
found in CSEs. We have used a method similar to that 
used in FG to determine the magnetic field strength par- 
allel to the line of sight. Like FG we are limiting ourselves 
to a Zeeman interpretation of the observed splitting, al- 
though a non-Zceman interpretation has been presented 
in Ncduloha & Watson (1990). 

2. Observations 

The observations were performed at the NRAOJj Very 
Long baseline Array (VLBA) in December 1998. We ob- 
served 4 late type stars (S Per, U Her, VY CMa and NML 
Cyg), the results presented here are the first results for the 
supergiant S Per, the source with strong water maser fea- 
tures and relatively simple structure. The stellar velocity 
Vlsr of S Per is -38.5 km s'^The beam width at 22.235 
GHz, the frequency of the 616 — 623 rotational transition 
of H2O, is « 0.7 X 0.3 mas. This allows us to resolve the 
different H2O maser features in the CSE. The data were 
correlated twice, once with modest (7.8 kHz = 0.1 km s^^) 
spectral resolution, which enabled us to generate all 4 po- 
larization combinations (RR, LL, RL and LR). The second 
correlator run was performed with high spectral resolution 
(1.95 kHz = 0.027 km s~^), necessary to detect the cir- 
cular polarization signature of the H2O Zeeman splitting, 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
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and therefore only contained the two polarization combi- 
nations RR and LL. The data produced by the first corre- 
lator run were used to accurately calibrate the R- and L- 
polarization. The calibration solutions were obtained on R 
and applied to L after we determined the relative R-L cor- 
rections, assuming the R- and L-polarization line profiles 
to be similar. The solutions were then applied to the data- 
set produced by the second correlator run, which we used 
to determine the circular polarization V for the separate 
maser features. The data analysis followed the method of 
Kemball, Diamond & Cotton (1995). 

3. Background 

The H2O (6 16 — 523) rotational transition consists of 6 hy- 
perfine components. Analyses of interstellar water masers 
indicate that all 6 hyperfine components contribute to the 
maser (Walker, 1984). The observed intensity ratios how- 
ever, frequently deviate from those obtained from molecu- 
lar transition probabilities (Moran et al., 1973; Genzel et 
al., 1979). Our analysis is performed using fitted hne ratios 
for the 3 strongest hyperfine components {F = 7 — 6,6 — 5 
and 5 — 4). The separation between these components is 
0.45 and 0.58 km s^"'^ respectively. The weakest compo- 
nents (i^ = 5 — 6, 5 — 5 and 6 — 6) are separated by more 
than 2.5 km s^^. They are not observed in the total power 
spectrum so we have not included them in our analysis. 
Our treatment follows closely the analysis performed in 
FG. Here we have added the possibility of multiple maser- 
ing hyperfine components. We assume the strongest hy- 
perfine component (i^ = 7 — 6) to be the dominant transi- 
tion. Because of the non-paramagnetic nature of the H2O 
molecule, the Zeeman splitting is extremely small. It is due 
to the interaction of the nuclear magnetic moment with 
the external B field. Thus the splitting is a factor of IQ-^ 
weaker than that for radicals like OH. In the weak field 
limit, the split energy AEz of a given energy level (J, F.I) 
is determined by: 



AEz 

with : 

aj 

aj 



-{ajgj + aigi} ■ ^J,NMF ■ Bca 



(1) 



{J{J + 1) + F{F + 1) - /(/ + 1)}/2F{F + 1), 
{F{F + 1) + /(/ + 1) - J{J + 1)}/2F{F + 1). 



This corresponds to a characteristic frequency shift of the 
order of Ai^z « 10^ Hz-[SGauss]- Where /(= 1) is the 
nuclear spin, fiN the nuclear magneton and Mp the mag- 
netic quantum number; gj = 5.585, and the jij-factors 
56 = 0.6565 and 35 = 0.6959 are from Kukolich (1969). 

Each hyperfine component will split into 3 groups of 
lines (cr+ , a~ and tt) , as seen in Fig. H for one of the hyper- 
fine components. The relative strengths of the transition 
probabilities have been determined by Deguchi & Watson 
(1986). For a magnetic field B parallel to the line of sight 
the Zeeman pattern consists of the two circular polarized a 
components only. The right- and left-handed (RR and LL) 
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Fig. 1. The Zeeman pattern of the i^ = 7 — 6 hyperfine 
component for an external field of 50 mG. 
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Fig. 2. Synthetic V-spectra for i^ = 7 — 6,6— 5 and 5 — 4 
calculated for an external field of 50 mG and a line width 
of Awl = 0.5 kms-i. 



spectra, corresponding to the a components will only 
be slightly shifted against each other (Awz ~ 10^"^ '° ^"^ 
times Avh)- As a result, the observed V-spectrum (RR- 
LL) will be a sine-shaped function, corresponding to the 
derivative I' of the total power spectrum. The ampli- 
tude of this function depends on the maser line width, 
the magnetic field strength, and on which hyperfine com- 
ponents actually contribute to the maser. By calculating 
synthetic V-spectra from the Zeeman pattern for different 
line widths, magnetic field strength and hyperfine combi- 
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nations we find the following relation for the percentage 
of circular polarization: 



Pv 



(K, 



*^niin) / ^n 



Av- 



F-F' 



B. 



Gauss 



/Ai^L. 



(2) 



Here V^ax and Vmin are the maximum and minimum 
of the synthetic V-spectrum fitted to the observations. 
At)L[kni s^ ] is the Gaussian line width of the total power 
spectrum, and /max is the peak flux of the maser fea- 
ture. B is the magnetic field component along the line 
of sight. The Ai?_i?'[-10^'^] coefficient depends on the 
masering hyperfine components. The Ap-pi coefficients 
have been determined by calculating Py from synthetic V- 
spectra, determined for a series of magnetic field strengths 
B, and for the different hyperfine components. For the 
F = 7—6,6 — 5 and 5 — 4 components individually we find 
Ap^F' = 16.22, 10.00 and 1.23 respectively. For a fitted 
combination of hyperfine components we find slightly dif- 
ferent values. An example of synthetic V-spectra for the 
three hyperfine lines is shown in Fig. g. 

However, due to the complex interactions between the 
various hyperfine components in the maser regime, de- 
viations from the V-spectrum proportionality are possi- 
ble. Detailed radiative transfer treatment, as performed by 
Nedoluha & Watson (1992; hereafter NW), for instance, 
resulted in Af-f' ~ 23.9. 

4. Results 

Fig. shows the total intensity map of the water maser 
features surrounding S Per. We are able to identify most 
of the maser features detected in earlier observations (Di- 
amond et al., 1987; Marvel, 1997). The positions are rel- 
ative to the brightest maser feature, for which we have 
managed to determine the circular polarization spectrum 
shown in Fig |j. This figure also shows a x^-?A to the sine- 
shape spectrum. The amplitude of the V-spectrum is only 
a small fraction (« 1%) of the total power, so we have only 
been able to determine the Zeeman splitting, and thus 
magnetic field strength, for the brightest maser feature. 
We do not detect circular polarization in any of the other 
bright maser features although, if it was present at the 
same absolute level as in the brightest feature, we would 
have detected it. This further confirms our detection, be- 
cause a scaled down version of the total power / would 
also have been detectable on the other strong features if 
calibration errors were significant. 

The total power spectrum indicates that one of the hy- 
perfine components clearly dominates, since the splitting 
of the hyperfine lines should otherwise have been observ- 
able. We have performed a fit to the total power spectrum 
to determine the maser line width and the best fitted ra- 
tio for the three strongest hyperfine transitions. Using this 
line ratio we have calculated the Af-f' coefficient as de- 
scribed above. In this case we find Af-f' = 15.54. Us- 
ing the fitted Gaussian line width (At;L = 0.44 ± 0.01 
km s-i) and Py = (9.9 ± 0.5) • 10"^ in Eq. |, we find 
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Fig. 3. Total intensity image of the H2O maser features 
around S Per. 



for the magnetic field strength along the line of sight 
B\\ =279±30mG. Ifonly theF = 7-6, 6-5 or 5-4 hy- 
perfine transition contributes the magnetic field should be 
scaled by 0.96, 1.55 or 12.63 respectively. As seen in Fig. H 
the V-spectrum is negative on the blue shifted side, there- 
fore the observed component B\\ is pointing away from 
us. This result is the first measured circular polarization 
in the circumstellar water maser region. 

5. Discussion 

As discussed above, the magnetic field strength was ob- 
tained by using a best fitted line ratio for the three main 
hyperfine components. A radiative transfer treatment for 
the polarized maser radiation of the 616 — 523 H2O ro- 
tational transition was performed by NW. Some difficul- 
ties exist matching their results to our observations of the 
total intensity spectrum, which shows an almost Gaus- 
sian shape. Their treatment also did not predict the anti- 
symmetric shape of the V-spectrum as shown in our ob- 
servations and those by FG. However, since the calibra- 
tion performed here and by FG assumes similar R- and 
L-polarization line profiles, the anti-symmetric shape is a 
necessary result of the treatment of the data. 

The observed V-spectrum however, is not directly pro- 
portional to /' either, as would have been expected in the 
simple model. The minimum and maximum of the sine- 
shaped function are not located at ±cr/v2, with a being 
the observed Gaussian line width. Instead they are found 
at Aj/ smaller by a factor of « 2.5 with respect to the 
line width fitted to the total power spectrum. Possibly, 
this narrowing of the V-spectrum can be attributed to 
the overlap of the multiple hyperfine components, as pre- 
dicted by the treatment and analysis of NW. The observed 
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Fig. 4. Total power (/) and circular polarization (V) spec- 
trum of the brightest H2O maser feature around S Per. 
The dashed line is the fit of the synthetic ^-spectrum 
to the observed spectrum. Also shown are the observed 
(dashed) and expected (solid) positions of the minimum 
and maximum of the V-spectrum. 

effect, however, still seems too large. Due to this our mag- 
netic field strength could be overestimated by at most a 
factor of 2. The interstellar water masers observed by FG 
did not show this narrowing. 

Because of the narrowing of the V-spectrum it is diffi- 
cult to address the saturation state of the maser. Elitzur 
(1998) showed that the observed circular polarization 
spectrum would be a good indication of the maser sat- 
uration state. If the ratio of V{i')/I'{i^) increased towards 
the wings of the line instead of showing a constant ratio 
the maser is thought to be unsaturated. This is however, 
almost the opposite effect of what we observe. Until the 
narrowing of the circular polarization is fully explained 
our observations are difficult to reconcile with a specific 
saturation state. 

We have also tried to determine the linear polarization 
of the maser features around S Per, but no indication of 
linear polarization has been found. This is consistent with 
the observations of Barvainis & Deguchi (1989). They ex- 
plain that the absence of linear polarization is probably 
due to the fact that masers are not very strongly sat- 
urated, and that infrared radiation does not contribute 
significantly to the pumping process. 

The magnetic field strength derived this way is within 
the range estimated by previous observations of SiO and 
OH masers. As noted before, polarization observations of 
SiO masers close to the central star reveal fields of 5-10 
G, assuming a standard Zeeman interpretation. OH maser 
observations of features around S Per indicate a field of 
slightly less than 1 mG (Masheder et al., 1999). Based on 



these values, the dependence of -B oc r~^ for a solar-type 
magnetic field is the most likely. For a dipole medium, the 
magnetic field is expected to vary with r^'^, which appears 
to be too steep to accurately describe the observations. 

In S Per, the H2O and OH masers are observed to ex- 
ist at similar projected distances (Masheder et al., 1999). 
This would disagree with the observed differences in the 
magnetic field strengths, except if the magnetic field could 
remain frozen in high density clumps. The magnetic field 
strength is then expected to vary with number density 
a.s B (X n'^, with 1/3 < k < 1/2 from theoretical pre- 
dictions (e.g Mouschovias 1987), where n is the number 
density. SiO masers are observed in high density clumps 
at 5 — 10 AU from the central star. H2O masers exist in 
similar clumps at distances of a few hundred AU, with the 
magnetic field lines frozen in the dense medium. Richards 
et al. (1999) show that the OH and H2O maser clumps 
avoid each other, although located at similar projected 
distances. They conclude that the density ratio between 
the H2O maser clumps and the OH in the surrounding 
medium only needs to be a factor of 50. However, mag- 
netic fields frozen into the maser clumps would require a 
density ratio of « 10** to explain the difference in field 
. This seems to indicate that actual coexistence between 
the OH and H2O masers is unlikely. 

In conclusion, although the exact infiuence of the hy- 
perfine interaction is not yet clear, we derive a magnetic 
field strength of B|| = 279 ± 30 mG. 
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